The Blazar Main Sequence 

A. Cavaliere, V. D'Elia 
Astrofisica, Dip. Fisica 2a Univ. Tor Vergata, Roma 1-00133 



Received 



accepted 



- 2 - 



ABSTRACT 

We propose a sequence (the Blazar main sequence, BMS) that hnks the 
two main components of the Blazar class, namely, the Flat Spectrum Radio 
Quasars and the BL Lacertae objects, and yields all their distinctive features 
in a correlated way. In this view, both type of sources are centered on a 
supermassive Kerr hole close to maximal spin and observed pole on. But the 
FSRQs are energized by accretion at rates m ~ 1 — 10, and are dominated by 
the disk components (thermal and electrodynamic jet-like component) which 
provide outputs in excess of L ~ 10^^ erg s"^. On the other hand, accretion 
levels m <^ 1 are enough to energize the BL Lacs; here the radiation is highly 
non-thermal and the power is partly contributed by the rotational energy of 
the central Kerr hole, with the latter and the disk together sustaining typical 
L ~ 10^^ erg s~^ for several Gyrs. If so, we expect the BL Lacs to show quite 
different evolutionary signatures from the FSRQs, and in particular, number 
counts close to the Euclidean shape, or flatter if the sources make a transition 
to the BL Lac from a FSRQ mode. In addition, for lower rh along the BMS we 
expect the large scale electric fields to be less screened out, and to accelerate 
to higher energies fewer particles radiating at higher frequencies; so in moving 
from FSRQs to BL Lacs these non-thermal radiations will peak at frequencies 
inversely correlated with the disk output. For the BL Lacs such dependence 
implies increased scatter when one tries a correlation with the total ouput. At 
its endpoint, the BMS suggests widespread objects that are radiatively silent, 
but still efficient in accelerating cosmic rays to ultra high energies. 

Subject headings: Quasars: general - quasars: evolution - galaxies: nuclei - 
galaxies: groups 
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1. Introduction 

The Blazars are commonly perceived as a class of Active Galactic Nuclei marked out 
from the rest of their kin by a number of features: strong, high frequency radio emission 
from compact cores, often with superluminal expanding components (see Jorstad et al. 
2001); powerful 7 rays extending to energies hu ~ 10^ MeV and beyond (see Mukherjee 
1999); rapid multi- wavelength variability (see Bottcher 2000); high and variable optical 
polarization (see Yuan et al. 2001). Comprehensive reviews and extensive references 
concerning these sources may be found in Urry & Padovani (1995), Padovani & Urry (2001, 
henceforth PUOl). 

The above features defining the Blazar class are widely explained in terms of the 
view originally proposed by Begelman, Blandford & Rees (1984). This holds the radiation 
from all these sources to be produced in a coUimated, relativistic jet of particles with bulk 
Lorentz factor F ~ 5 — 20; when observed at small angles of order F"^ the jet produces the 
"blazing" effects that provide the class denomination. The intrinsic luminosity emitted by 
the jet over the entire solid angle, and total jet kinetic power are reviewed by Ghisellini 
(1999a), whose definitions and beaming factors we will adopt in the following. 

Within the Blazar class two main subclasses are usually identified, namely, the BL Lac 
Objects and the Flat Spectrum Radioloud Quasars (FSRQs). The former differs from the 
latter subclass on several accounts, that we shall discuss in this paper: a) lack of emission 
lines with EW ^ 5 A, and lack of blue-UV bumps; b) total outputs L ^ 10^^ erg s~^, 
compared with FSRQs often exceeding this limit and in some cases approaching 10^^ erg 

(see specifically Maraschi 2001); c) spectral energy distributions peaking at frequencies 
from optical to X-ray and in the range of 10 GcV, compared with the FSRQ peaks at 
far IR to optical frequencies and around 10~^ GeV; d) no signs of cosmological evolution, 
compared with the strong evolution clearly displayed by the FSRQs in common with the 
other quasars. 

These distinctive features stand out of the many selections entangling this field, and 
even though intermediate objects are currently found (Sambruna et al. 2000, Perlman et al. 
2001), the FSRQ/BL Lac strong bimodality calls for explanation; this is the scope of the 
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present work, that expands some preliminary results reported by D'Elia & Cavaliere (2001). 

We will base on the accreting black hole (BH) paradigm (see Begelman, Blandford 
& Rees 1984), in the framework of a 4-parameter classification similar to that discussed 
by Blandford (1990, 1993); this comprises the BH mass M, its angular momentum J, 
the accretion rate M, and the viewing angle 9 from the axis of the jet (if any). In the 
scheme wc contemplate, radio quiet AGNs correspond to slowly rotating BHs with low 
(Seyfert galaxies) or full (quasars) M compared to the Eddington rates; the quasars require 
large BH masses M of order 10^ Mq to attain proportionally high Eddington luminosities 
Le ^ lO^^Mg erg s-^ 

Radio loudness constitutes a debated issue. Radio loud AGNs are often observed to 
reside in very bright galaxies. On the one hand, these sources have been inferred to correlate 
with large BH masses; for instance, Laor (2000) stresses the case for radio loudness being a 
consequence of large M. On the other hand, the majority of the optically selected quasars 
including very bright ones are radio quiet (see for recent data Goldschmidt at al. 1999, 
Hamilton et al. 2001); in addition, Scarpa & Urry 2001 find that the hosts of quasars strong 
in the radio are no different from giant ellipticals in general, including hosts of radio quiet 
quasars. Finally, Ho (2001) from a comprehensive database finds no simple relation between 
radio loudness and M. Similar data from more limited observations spurred the view (see 
Blandford 1993) that the other BH parameter, namely the spin, should be involved as a 
necessary, intrinsic condition for directing and launching the jets associated with radio 
emissions. In the present paper we focus on this view. 

In fact, we will consider the Blazars in the context of supermassive holes with 
M ~ IO^Mq close to the maximal observed values for MDOs (see Gebhardt et al. 2000, 
Ferrarese & Merritt 2000), with high angular momentum J ~ Jmax close to the maximal 
Jmax = GM'^/c, observed nearly pole on within angles 9 ^ F"^ (see Sambruna, Maraschi & 
Urry 1996). We shall concentrate on discussing how all Blazar properties can be arranged 
into a Blazar main sequence (BMS) in terms of levels of the accretion rate m = M /Le- 
Levels m ~ ri~^ ~ 1 — 10 (the radiative efficiency being up to about 0.4 around a Kerr hole) 
mark the FSRQs, while values m ^ 10~^ mark the BL Lacs. 
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Our plan is as follows. In §2 we will stress how the optical properties of the Blazars 
can be understood in terms of different levels of rh. In §3 we show that high rates m ~ 1 
can produce the huge thermal and non-thermal outputs featured by many FSRQs, while 
conditions where m <^ 1 lead to the moderate, but mainly non-thermal outputs of the BL 
Lacs. In §4 we will argue that for lower values of m along the EMS the typical particle and 
photon energies are expected to increase; a secondary physical parameter relevant at low 
m causes considerable variance within the BL Lac subclass, adding to the effects from the 
orientation. In §5 we show how in the same framework we understand the widely different 
evolutionary properties observed in FSRQs and in BL Lacs. Finally, in §6 we summarize 
and discuss our conclusions. 



2. Optical Properties 

In the standard view (see Frank, King & Raine 1992; Peterson 1997) all thermal 
outputs are produced in, or around the accretion disk, with only weak anisotropics and 
with power given by Lth ^ ijTfiLE. 

Here we take up this view, and just stress how from the broad optical emission lines 
and the blue-UV bump in particular wc derive a first indication toward a BMS in terms 
of levels of rh. The bump is produced by power reprocessed in a superposition of black 
body emissions from the inner rings of the accretion disk, or in a hot corona (see Sun & 
Malkan 1989, Siemoginowska et al. 1995); the lines are produced by the high UV-soft X-ray 
continuum reprocessed by distant gas "clouds" (see Netzer 1990). In both cases the density 
of the reprocessing material is expected to scale up with increasing rh. 

The FSRQs share with other quasars conspicuous bumps and prominent broad emission 
lines; the strength of these features are consistent with high accretion rates m ~ 1. 

On the other hand, the weakness or absence of such features in the BL Lac spectra 
can be understood in terms of m ^ 10~^. In fact, at such levels of rh not only a weaker 
ionizing/exciting continuuum is emitted by the disk, but also lower gas densities are 
expected in the disk and its environment out to pc and larger scales; these conditions 
concur to account for the weak or absent emission lines. 
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The other effect that in the BL Lacs acts to swamp any residual bump and to further 
reduce the EW of intrinsically weak emission lines is constituted by the blazing continuum 
from the jet when observed at a small angle, see Vagnetti, Giallongo & Cavaliere (1991); 
Georganopoulos (2000). 

We recall that low values of m related to radiosources had been proposed by Bcgclman 
et al. (1984). Specifically, the different power levels and optical spectra of the FRIl and 
FRI radiosources have been previously related to high and low values of m, respectively, by 
Baum, Zirbel, & O'Dea 1995; this point will be further discussed in §6. 

3. The Blazar Power 

As to the jets, we assume them to be powered by variants of the mechanism originally 
proposed by Blandford & Znajek (1977, BZ) for direct extraction of rotational energy from 
a Kerr hole via the Poynting-like fiux associated with the surrounding magnetosphere. The 
BZ power scales as 

Lk ^Blrl, (1) 
in terms of the magnetic field threading the hole horizon at Vh- 

Variants are necessary in view of the limitations to recently stressed by Modersky 
& Sikora (1996); Ghosh & Abramowicz (1997); Livio, Ogilvie & Pringle (1999). Such 
variants involve either high strengths of Bh, or the similar electrodynamic output of the 
disk from a larger radius ra. 

In the first variant, very strong Bh as advocated by Meier (1999) arc required to 
account for the huge outputs of some FSRQs. Fields up to -B^/Stt ~ pc^ in the plunging 
orbit region have been argued by Krolik (1999); Armitage, Reynolds & Chiang (2000) and 
Paczynski (2000) discuss why such field values are unlikely in a thin disk. In thick disks the 
status of such enhanced fields is still unsettled. 

The second variant takes up from recent discussions (Ghosh & Abramowicz 1997; Livio, 
Ogilvie & Pringle 1999) that have stressed the continuity of Bh with the field Ba at the 
inner rim of a standard a-disk; in turn, i?^ is bounded by the maximum pressure, following 
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B^/An < Pmax- Depending on m, the inner disk is dominated by gas pressure (GPD) or 
by radiation pressure (RPD); the pressure scales as Pmax ^ {aM,^)~^/'^^m^_ll in the first, 
and as as Pmax oc (aMg)"^ in the second regime. Correspondingly, in the first regime the 
maximum power extractable from the hole grows with m, and in the jet frame the emission 
reads (see Modersky & Sikora 1996; Ghosh & Abramowicz 1997) 

Lk = lQ''Ml''''m1l{J/Jmaxf erg s"^ {GPD) (2) 

in terms of the Kerr hole mass M — Mg 10^ Mq; but for m ^ 10~^ the power Lk saturates 
to 

Lk = 2 10^^ Mg {J/Jmaxf erg s"^ {RPD) . (3) 

This is because the radius Tc defining the RPD region grows with m^^/^^ (Novikov & Thorne 
1973), and it turns out to exceed the radius r^s of the last stable orbit when m ^ 10~^ 
holds; in the latter conditions an inner RPD region actually exists. 

On the other hand, the electrodynamic power L^ contributed by the disk and itself 
strongly anysotropic scales as 

La oc Bl rl ; (4) 

here the effective disk radius Vd is bounded at several times Vh either by the value of 
Tf. or by the radial decline in the disk of the efficiency f] ~ G M M /2r . Wc recall that 
the ratio L^/Lk had originally been estimated by Blandford & Znajek 1977 at values 
^ 5; recent reappraisals (e.g., by Livio, Ogilvie & Pringle 1999) have evaluated it at 
La/LK ~ {Ba/Bn)^ (r./r,)^/^ 

We draw the following implications concerning the BMS. Many BL Lac jets can live on 
accretion rates m ~ 10~^, since the power levels Lk given by eq. (3) are adequate even 
including the kinetic energy in the jets (Celotti, Padovani & Ghisellini 1997, Maraschi 
2001). A larger but comparable contribution comes from the disk, so that Lk ^ L^^ Lth 
holds. The disk is coupled to the rotating hole both dynamically and magnetically (see 
Livio, Ogilvie & Pringle 1999) through the magnetic field fines threading the hole and 
tethered in the disk. 

On the other hand, many FSRQs feature total outputs L > IQi^^ erg s~^, with specific 
sources exceeding lO^^erg s~^, see Tavecchio et al. (2000), and Maraschi (2001), who include 
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a high proton contribution to the jet energy and adopt average values of F ^ 10 (see Sikora 
2001). Such luminosities require not only a saturated (see eq. 3), but also a dominant 
component up to Ld ~ 10^ Lk from the disk corresponding to Bd a few times larger than 
Bh and ~ Sr^; this implies an extended region dominated by radiation pressure, hence 
conditions where m ~ 1, consistent with the requirements from the optical properties 
discussed in §2. 

We note that in FSRQ conditions both Ld and the thermal emission Lth are fed by the 
gravitational power supply of order G M M jlvd- The share between and Lt^ depends 
on the power spectrum of the magnetic field inhomogeneities, that is, on the power share 
between two scales: the large scale, coherent vs. the small scale, turbulent component. A 
dominant share of the former, such as to yield a large Ld and a large associated transfer of 
angular momentum outward, may require modifications of the standard a-disk model as 
discussed by Salvati (1997). 

The hole output L^, even when it constitutes a minor component to the total, is likely 
to be important to provide on the jet axis a "high-velocity spine" instrumental for the 
outward jet propagation, see Livio (2000), Chiaberge et al. (2000). 

4. The Blazar SEDs 

Another feature of the Blazars is constituted by their extended spectral energy 
distribution (SED). This shows a first peak commonly interpreted as synchrotron radiation 
(see Urry & Padovani 1995). In the FSRQ spectra the peak frequencies v^eak he mainly in 
the range from the far IR to the optical band; in the BL Lacs the peaks are generally shifted 
to higher frequencies, and in some objects up to the hard X rays (Costamante et al. 2001). 

An additional high energy component to the radiation (see PUOl) is observed into the 
10 GeV range for the FSRQs, and out to several TeVs for the BL Lacs Q. This is likely to 

^Note that the 7-ray spectra at the sources may extend to higher energies than observed 
at Earth after the intervening absorption due to the IR background. If so, the following 
argument concerning coherent acceleration would be reinforced 
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be inverse Compton emission produced by substantial numbers of GeV or of 10^ — 10^ GeV 
electrons, respectively, scattering off external or synchrotron seed photons as discussed by 
Ghisellini (1999a). 

The apparent inverse scaling of Vpeak with the radiative L in moving from FSRQs to 
(quiescent) BL Lacs has been pointed out by Fossati et al. (1998), and interpreted as 
a cooling sequence due to the radiative energy losses being faster in the more powerful 
sources, up to controlling the maximal electron energies. Concerning the acceleration, one 
view considers weak electric fields over relatively large distances; these occur in the internal 
shock scenario where, however, the typical energies produced are of order £^ ~ 2 F < 20 
GeV and fall short of the top energies required in some BL Lacs as discussed by Ghisellini 
(1999b). 

So we are lead to consider as an addition or as an alternative particularly relevant to 
BL Lacs, the electron acceleration due to higher electric fields acting over shorter distances. 
We shall see that this leads to an acceleration sequence also linked with the BMS; this is 
because in the BL Lacs lower levels of m allow for higher electron energies and so produce 
a SEDs peaked at higher frequencies. Meanwhile, a 2"*^ intrinsic parameter produces the 
specific scatter of I'peak within this subclass. 

Such high E fields naturally arise around the BZ magnetosphere in association with 
considerable energy transport via "Pointing fiux" along the jet, with the magnetic field 
decreasing outwards rather slowly, like B = Bd (rms/r)^'^^^^ with p < I; typically p = 1/4, 
in a similar manner to Blandford & Payne (1982). The force-free condition E • B — 
governing the BZ magnetosphere has to break down at the flow boundaries, and to give 
way to electric fields E ^ B; this is likely to occur time-dependently and inhomogeneously 
within sheets or filaments located at average distances r ~ 10^^ -^ 10-^^ cm. But such 
E fields arc clcctrodynamically screened out beyond distances exceeding a few times 
d — c/up oc (7/n)^/^, where Up — (47re^n/m7)^/^ is the plasma frequency of the screening 
particles with Lorentz factor 7 and density n. 

Lower bounds to the densities may be estimated from the powers emitted in the form 
of synchrotron or inverse Compton radiation, that is, L ^ ^'^U E? n where U is the total 
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energy density in quiescent conditions. We obtain n ~ 1 cm~^ for BL Lacs with L ~ 10^^ 
erg s~^, and n ~ 10^ cm~^ for FSRQs with L ~ 10^^ erg s~^. Correspondingly, we expect 
many electrons to attain high energies, up to 

-i^ax ^eBd/nieC^^ 10' ^4 dio rr/-'' , (5) 

where we have used d — 10^° dio cm, and r — 10^^ ri7 cm; we have also expressed as 
Bd = 10^ B4 G the field values at r^s ~ 1-5 10^"^ Mg cm, which decrease into the emitting 
region like B{r) = B^ {rms/f'Y''^^ as recalled above. 

In eq. (5) there is room for variations of B and n; the former may be down by a 
factor 10 or the latter up by factors 10^ (with M again set at lO^M©), and eq. (5) still 
yields 7 ~ 10^. This is enough to produce inverse Compton photons up to ~ 10 TeV, the 
highest energy photons observed to now from the BL Lacs. In FSRQs densities 10^ larger 
give electron energies some 30 times smaller, consistently with the parameters derived by 
Bottcher & Dermer 2001 from advanced spectral modeling. Given ''jmax^ we expect less 
energetic but more numerous electrons to arise naturally when the particles escape before 
full acceleration or E < B holds, and when the inhomogeneities of B are considered; the 
latter varies outwards following B oc r~^'^^ as said, and this condition by itself produces an 
energy range by factors close to 10^ within a distance range r ~ 10^^^^ cm. 

The electron energies given by eq. (5) are proportional to d oc and so scale hke 

7„.„. oc (7^C/i?3)V2 (6) 

This results in 7,^,01; oc L~^/^ if we adopt the empirical relation 7^C/ ~ const reported by 
Ghisellini (1999a) for the FSRQs. On the other hand, the same author finds that the 
7 values for the BL Lacs exceed this relation, but may be included on modifying it to 
read 7^C/ ~ [/~°-^; from this we obtain '-)rnax L~^/^. Thus we expect the peaks of the 
non-thermal emissions, and more clearly the syncrotron Vpeak oc 7max5 anticorrelate 
strongly with the disk luminosity following Vpeak L~^/^ or steeper, when comparing the 
FSRQ with the BL Lac subclass. 

In closer detail, the emitting particles are likely to constitute the high energy subset of 
the total particle population; but similar results obtain if all electrons are described as in 
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Ghisellini (1999a) by a power-law energy distribution ni^'j) oc " from ^rnin ~ 

10=^; with 

a ~ 2.5 we obtain Up^ak oc 7^ax 1'"^''^, while a steeper relation t'peafc oc holds with 
a ~ 2. Concerning the full SED wc refer to the complete modeling of the emitting regions 
by Bottcher & Dermer 2001, who treat quantitatively the spectral shapes radiated by 
synchrotron and inverse Compton emissions in the jets, considering also the cooling limit to 
the electron energies provided by the high photon densities in FSRQs; they evaluate specific 
parameters agreeing with the orders of magnitude we estimate. 

The other outcome from our view is that within the BL Lac subclass we also expect 
the second parameter Lk/L^ to affect the scaling laws provided by eq. (6). This is because 
the total radiative luminosity of the BL Lacs L = + Lk comprises a disk and a Kerr 
hole component with comparable powers. Of these, only ~ G M M /2dms is directly 
related to m and to n while Lk, when considerable, tends to be independent of them as 
shown by eq. (3). Thus at a given (hence at given M and n) higher total luminosities 
obtain in BL Lacs where Lk is important. In other words, Vpeak scales inversely with the 
component La only, and it is bound to show scatter specific to the BL Lacs when one tries 
to correlate it with the total luminosity. As a consequence, we do not expect the scaling to 
extend unblurred from the so called LBLs to the HBLs (low and high energy peaked BL 
Lacs, respectively, see Padovani & Giommi 1995). In fact, considerable scatter of Vpeak a-t 
given L is increasingly found within the BL Lac subclass, see Giommi et al. (2001). 

On the other hand, for the FSRQs L ^ holds, that is, Lk is irrelevant as to the 
total power, and so is its contribution to the scatter. 

The specific BL Lac scatter discussed here adds to other components common to all 
Blazars due to different values of -B, M, and viewing angles (see Sambruna, Maraschi & 
Urry 1996, Georganopoulos 2000), and to variations of Up^ak and L arising during flares. 

To this point we have shown how all spectral properties of FSRQs and BL Lacs are 
expected to consistently vary along the BMS. In the next section we discuss why we expect 
the Blazar evolutions to be also related to different levels of m. 
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5. The BlciZcir Evolutions 

Strong cosmological evolution is closely shared by the FSRQs with the rest of the 
quasars (see Wall & Jackson 1997, Goldschmidt et al. 1999), and shows up, e.g., in their 

steep number counts; the BL Lacs, instead, show no signs of a similar behavior (Giommi, 
Menna & Padovani 1999; Padovani 2001). 

Here we argue that such an evolutionary peculiarity is related to the following two 
independent circumstances: i) At levels m ^ 10^^ all BL Lacs shine for several Gyrs with 
only a mild decrease in luminosity (a very slow "luminosity evolution", LE), and this 
implies counts close to Euclidean, ii) A genetic link is conceivable from FSRQs to BL Lacs, 
that is, some of the former may switch from a short-lived regime of high rh to the long-lived 
regime of low m; so the BL Lacs undergo also some negative "density evolution" (negative 
DE), and their counts are further flattened. Since point i) has been discussed by Cavaliere 
& Malquori 1999, we just recall their results below, and concentrate on point ii). 

The evolution of the bright optically selected (and mainly radio quiet) quasars is 
strong, see Boyle et al. (2000), and so is in X-rays, see Delia Ceca et al. (1994); Miyaji, 
Hasinger & Schmidt (2000). For the bright sources, this behavior is widely traced back 
to the exhaustion in the host galaxy of the circumnuclear gas stockpile usable for further 
accretion onto the central black hole. Such an exhaustion is caused by previous accretion 
episodes and by ongoing star formation (Cattaneo, Haehnelt & Rees 1999; Kauffmann & 
Haehnelt 2000; CavaUere & Vittorini 2000, CVOO); the ensuing decrease of the average 
accretion rate rh is fast, and causes the average bolometric luminosities to scale down 
rapidly for z < 2.5 (a strong LE). The radio loud quasars share this evolutionary behavior, 
constituting a fraction (not yet quantitatively understood, see §1, and the discussions by 
Moderski, Sikora & Lasota 1998 and by Livio 2000) of order 10% of the total; the FSRQs 
in particular are widely held to be the same sources seen pole on. 

But as the rates m decrease we also expect a considerable fraction of the underlying 
hole-disk systems to eventually switch from being overwhelmingly fueled by accretion to 
being partly fed by the hole rotational supply Ex] the latter had been stockpiled during 
the phases of accretion of coherent angular momentum J along with the mass (Bardeen 
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1970; Moderski, Sikora & Lasota 1998). In other words, we expect some radio loud sources 
to switch from the FSRQ to the BL Lac mode. The moderate outputs typical of this mode 
can be sustained for several Gyrs by the dynamically and magnetically coupled system 
constituted by Kerr hole and inner disk; so the decrease of BL Lac luminosities is expected 
to be slow (Cavaliere & Malquori 1999), with time scales in the range estimated to a first 
approximation from ~ Ek/Lk ~ 5 -j- 10 Gyr, rather uniform within the population. 

In closer detail, CVOO (see also refs. therein) trace the accretion episodes feeding the 
quasars back to interactions of the host galaxies with their companions in a group. These 
events destabilize the host gas, and trigger accretion over times of order r ~ a few 10~^ 
Gyr that comprise the duration of a close interaction and a few galactic dynamical times 
for the inflow to develop and subside; some 3 — 5 repetitions are expected per host galaxy 
since z ~ 2.5. The frequency of such interaction episodes diminishes in time, and gives rise 
to a positive, weak DE of the quasar population over scales td ~ several Gyrs. Meanwhile, 
the efficiency of such episodes drops as anticipated above, due to the exhaustion of the host 
gas; the latter is halved on the scale tl ~ 3 Gyr. The result is that the average luminosities 
drop over the same scale, in other words, a strong LE is produced in the quasar population. 

On this basis CVOO compute the luminosity function (LF) for the optically selected 
objects at 2; < 2.5 to read 

Here Nq(z) is the space density of groups of galaxies, t~^{z) oc (1 + z)^^^ is the average 
interaction frequency for a host galaxy in a group, and / = L/Lb{z) is normalized using 
the break luminosity that scales as Li,{z) oc (1 + z)^ in the critical universe; we use this 
simple scaling for numerical evaluations, since CVOO show that the observables do not vary 
much in the Qo + — ^ cosmology. For the FSRQs we will use here eq. (7) with the 
proportionality factor of order 5 10~^. 

Our main point here is that the episodes of powerful FSRQ activity feeding on high 
rh will last some 10~^ Gyr and die out over similar scales; the sources may be no longer 
reactivated after a "last interaction". This leaves behind a rapidly spinning hole (and a 
relatively quiescent though bright galaxy, see Scarpa 2001) in about 1/2 of the objects. 
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Thus the DE timescale td for the deaths of the bright FSRQs also constitutes the timescale 
for BL Lac births; we evaluate it from eq. (7) integrated over L > Lf„ to obtain 

^D-J^ln[j^ N(l,t)dl]^-ln[NG(t)/rr(t)], (8) 

with the second equality holding because only the time-dependent amplitude in eq. (7) is 
relevant here. We obtain values of td between 7 and 5 Gyr, depending on the form for 
small z of Ng{z) which may scale as (1 + z) on assuming a Press & Schechter (1974) mass 
distribution, or as (1 + z^ following Lacey & Cole (1993). 

One sign of evolution is provided by the integrated source counts, which may be 

evaluated from the relation (see Cavaliere & Maccacaro 1990) 

N{> S) oc 5-3/2[l - CiSo/Sy/^ + 0(3-')] . (9) 

This simple relation is limited to high-intermediate fluxes. On the other hand, it has the 
advantage of exposing how the evolutionary properties of the two Blazar subclasses depend 
on the two population time scales td and tl (the latter being long for the BL Lacs, and 
short for the FSRQs); these appear in the coefficient 

C = 3Do (/2)[2(1 + a) - 1/HoTd - (/? - 1)/ HoTl]/4 Rh{1'/') . (10) 

For the BL Lacs counted in the radio band a relevant value is r^, ~ 10 Gyr (marking the 
slow LE of the BL Lacs). Other quantities involved are: the flat spectral index a — —0.3 in 
the GHz range; the slope /3 2.5 of the radio LF; the normalized moments (Z") of the LF; 
the typical distance to high flux BL Lacs Dq = [Lq/ At^SqY^"^ Q.OSRh in Hubble units. 

On using only the slow LE on the scale tl 8 Gyr, eq. (10) with the Hubble constant 
Hq — 65 km/s Mpc yields C ~ — 0.1; then from eq. (9) the result is A'"(> F) oc F"^-^ as in 
Cavaliere & Malquori (1999). 

But BL Lac births imply also negative DE to occur, described by the scale td ~ — 6 
Gyr evaluated below eq. (8). When also this scale is introduced in eq. (10), the values of 
C increases to C Ri 0.1. Now the result is N{> S) oc S~^'^ at high-intermediate fluxes, 
consistent with the data by Giommi et al. (1999), see flg. 1; at 5 GHz similarly flat counts 
are both expected (see Jackson & Wall 2001) and observed (see the prehminary data by 
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Padovani 2001), while the FSRQ counts observed are quite steep (see Wall & Jackson 
1997, Jackson & Wall 2001). We add that no evidence of evolution has been observed by 
Caccianiga et al. (2001) from the values of Vg/K in a large sample of BL Lacs from the 
REX survey. 

The count slopes we evaluate for the BL Lacs are in sharp contrast with the result 
from a similar evaluation for the FSRQs. The parameters appropriate for them arc recalled 
above and read: td — + & Gyr, tl — 3 Gyr, while Dq fa 0.3 Rh holds and again P 2.5. 
Eq. (10) then gives C —8, and in turn eq. (9) yields radio counts A^(> S) steep well 
above S~^-^ and consistent with the data, see fig. 2. 

We thus see that also the evolutionary properties vary systematically along the BMS. 
Actually, the faint end of the LP is expected to be flattened by the varying beaming effects 
associated with a distribution of P values (Urry & Padovani 1995); this holds both for 
the FSRQs and the BL Lacs, and will contribute to the flattening of the counts at their 
respective faint ends, but it does not affect the intermediate range where sharply different 
slopes stand out. 

6. Discussion and Conclusions 

We have argued that Blazars constitute basically similar sources comprising a 
supermassive, nearly maximally spinning Kerr hole and its accretion disk; but large, 
systematic and correlated changes are bound to arise in all observational properties when 
the accretion rate levels vary from m ~ 1 that mark the FSRQs to m <C 1 that mark the 
BL Lacs. Such systematic changes basically arise because lower values of rh imply depleted 
densities in and around the disk; by the same token, these cause both higher energies of the 
emitting particles with less reprocessing of the emitted photons, and lower total luminosities 
over longer time scales leading to less evolution. 

Our view (in short, BMS for Blazar main sequence) outlined here lead us to specifically 
predict the following correlated features, also summarized in Table 1. 

a) Optical properties: in §2 we recall that high values of rh yield strong thermal 
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features, namely, emission lines and big blue bump. On the other hand, we stress that low 
values of rh imply intrinsic weakness or drowning of all these features, especially in the 
presence of a strong beamed continuum. 

b) Luminosities: in §3 we argue - in view of the limitations to the power directly 
extractable from a rotating hole - that disks with high values of rh are required to produce 
outputs as strong as those observed in FSRQs. In BL Lacs, on the other hand, low values 
of rh sustain intermediate or low disk luminosities, then comparable to those directly 
extractable from a Kerr hole after eq. (3). This yields a limiting power of about 10^^ erg 
s~^; it is intrinsic to the present view that the top BL Lac outputs should not considerably 
exceed such values, as in fact observed (see Maraschi 2001). 

c) SEDs: also the main features of the Blazar spectra can be predicted in this picture 
as shown in §4. Comparing BL Lacs with FSRQs, in the former sources we expect lower 
particle densities in the acceleration region, and less effectively screened electric fields. 
These produce higher electron energies up to some 10^ GeV, allowed on the other hand by 
the lower photon densities; the ensuing synchrotron emission will peak at higher frequencies. 
Prom b) and c) we derive the inverse correlation Upeak LJ^ or somewhat steeper, see §4. 
But in BL Lacs specific and considerable scatter is to arise when one tries to correlate Vpeak 
with the total luminosity L, because this is made up by two comparable components, Lk 
and Ld- 

d) Evolution: comparing FSRQs with BL Lacs, much weaker evolution is expected for 
the latter, as argued in §5. This is because their output is lower and is partly extracted 
from the Kerr hole rotational energy, in conditions of low and long lasting accretion rates; 
so these sources are are long lived, with a slow population evolution of the LE type over 
scales tl of several Gyrs. Using only this scale to predict the source counts from the 
relations given in §5, the result is A^(> F) oc F~^-^ at fluxes F > 10^^ Jy. In addition, 
we expect a genetic link to also occur; some FSRQs may make a transition to the BL Lac 
mode after a "last interaction" of their host with companion galaxies. Such BL Lac births 
are equivalent to a negative, though weak DE on a scale td ~ —6 Gyr; this competes with 
the slow LE to yield counts flatter yet, namely A^(> F) oc F^^'^ at intermediate fluxes, as 
shown fig. 1; equivalently, we expect values of K/K ~ 0.5, see Giommi, Menna & Padovani 
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1999. This scenario implies transitional objects to occur from FSRQs to BL Lacs, with 
m ^ 10~^ and shining with intermediate luminosities for a few 10^^ Gyr; in the optical 
band these sources will still show thermal features (emission lines and bumps), while their 
non-thermal SEDs will peak at higher Upeak compared with canonical FSRQs. On the other 
hand, similar relations computed with the shorter time scales appropriate for the FSRQs 
yield much steeper counts, as shown in fig. 2. 

We note that if the radiosources FRI and II constitute the parent populations of the 
BL Lacs and the FSRQs, respectively, then a genetic link between BL Lac births and FSRQ 
deaths implies a similar link between the related radiosurces. Such a link constitutes matter 
of current debate. We recall from §2 that power levels and optical spectra of the FRll and 
FRI radiosources have been related to high and low values of m, respectively, by Baum, 
Zirbel & O'Dea 1995. We stress that overlapping properties and truly transitional objects 
are being found (see Bicknell 1995, and PUOl with references therein); in addition, similar 
evolutionary behaviors are observed for the radiosources and for the Blazars related to 
similar values of m, that is to say strong evolution for the FRlI but weak or even negative 
evolution for the FRI radiosources (Jackson & Wall 2001). So there is scope for discussing 
"grand unification" schemes of radiosources and Blazars, as pointed out by PUOl. 

We also note that our BMS agrees well with the parallel work by Bottcher & Dermer 
2001; this is based on detailed modeling of the Blazar emitting region, and concerns the 
full SEDs evolving from FSRQs to BL Lacs as the accretion rates decrease and the black 
hole environments are depleted in gas, dust and reprocessed photons. In fact, the two 
papers turn out to be complementary, and concur toward the notion of a sequence or even 
a genetic link between the two subclasses. 

Finally, if we carry the above BMS to its extreme, that is, to residual m ^ 10~^, we 
are led to expect objects with very faint, if any, e.m. emission (further suppressed in the 
ADAF regime, see Di Matteo et al. 2000), along with nearly unscreened electric fields; these 
objects would be akin to the cosmic ray accelerators discussed by Boldt & Ghosh (1999). 

Under the widespread if arguable assumption that in such conditions a magnetic field 
B ~ 10^ — 10^ G can still be held around the hole with the slow upward decrease recalled in 
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§4, particles including protons would be accelerated over long distances, at radii that must 
exceed the force-free region's (acceleration at small radii has been shown to be untenable 
by Levinson 2000). These particles - in fact, cosmic rays with ultra high energies - are 
predicted by the analogous of eq. (5) but in conditions of unscreened fields, that is, large d 
due to very low n (and to consistently large 7). Thus we recover the long recognized energy 
limit given by 

Sma. -^eBar^, {r^s/rf/p ~ W B^mI'\^^'^ eV . (11) 

Tens of these accelerators could lie within 50 Mpc. As noted by Boldt & Ghosh (1999), 
they would evade in the simplest way the GZK cutoff at UHE; as to accounting for the 
observed particle flux, they need to produce only L ~ 10^^ erg s~^ mostly in energetic 
particles. Intergalactic magnetic fields of nG strength would blur to near isotropy the 
geometrical memory of the actual sources for most such UHECRs. 
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Table 1 - The Blazar Main Sequence summcirized 





FSRQs 


BL Lacs 


— > CR accelerators 


key parameter 


m ~ 1 


m ~ 10-2 


m ^ 10-^ 


optical features 


em. lines, bump 


~ no lines and bump 


none 


integrated power 


L ~ 10^^ erg s-i 


L < 10^5 erg s-i 


L < 10^2 erg 


Kerr hole vs. disk 


Lk < Ld 


Lk ^ La 


Lr > Ld 


top energies 


hu ~ 10 GeV 


hv ~ 10 TeV 


^max ~ 10^^ eV 


evolution 


strong 


weak if any 


negligible 



Figure Captions 



Fig. 1. The BL Lac counts evaluated from eqs. (8) and (9) using td — —7 Gyr 
(dashed) and td = — 5 Gyr (soUd). The dotted hne represents the Euchdean counts. The 
data points concern HBLs and are taken from Giommi et al. 1999, see their fig. 5. 

Fig. 2. The FSRQ counts evaluated from eqs. (8) and (9) using td — + 6 Gyr and 
tl = 3 Gyr (solid line). The open squares represent FSRQ counts from Wall & Jackson 
1997. 



-25- 




-26- 



~i 1 1 1 — I — I — r-r-i 1 1 1 — I — I I I I I 1 1 1 — I — I — I I i_ 



10 - 



1 r 



^Q-3 _ 



CM 

b) 

CD , 
CO 



\ 



□ 



_ I I I I I I I I I I I I I I I I I I I I I I I I I 

Radio flux (5 GHz) [Jansky] 



